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1. Introduction 


Entropy generation minimization or thermodynamic optimiza- 
tion is not an old technique. It is also components of exergy analysis 
[1]. The method of thermodynamic optimization or entropy gen- 
eration minimization (EGM) is a field of activity at the interface 
between heat transfer, engineering thermodynamics, and fluid 
mechanics. The entropy generation minimization technique is 
based on the simultaneous application of the first law and the 
second law in analysis and design. Entropy can be used to distin- 
guish between reversible and irreversible processes. It analysis can 
quantify the level of energy quality. In this paper we review the 
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fundamentals of the method, its current status, and a few examples 
on natural convection studies. 

The destroyed exergy is proportional to the generated entropy. 
Exergy is always destroyed, partially or totally: this is the state- 
ment made by the second law of thermodynamics. The destroyed 
exergy, or the generated entropy is responsible for the less-than- 
theoretical thermodynamic efficiency of the system. 

First law efficiency of thermodynamics in a heat transfer engi- 
neering system is very much restricted. Thus, calculations using the 
second law of thermodynamics, which is related to entropy gener- 
ation, are more reliable than first law based calculations. In almost 
all thermal systems, second law-based efficiency can be defined in 
terms of the ratio of actual thermal efficiency under the same condi- 
tions. Therefore, the second law of thermodynamics can be applied 
to investigate the irreversibility in terms of the entropy genera- 
tion rate. Determination of entropy generation is also important to 
enhance system performance, because entropy generation is the 
measure of destruction of available work of the system [2-4]. 
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Nomenclature 

Be Bejan number 

C concentration 

FTI fluid friction irreversibility 
HTI heat transfer irreversibility 
k thermal conductivity 

Ns total entropy generation 

S entropy 

T temperature 

u,v,w velocities 

x,y,z coordinates 

Se volumetric rate of entropy generation 
Greek letters 

H dynamic viscosity 

(o irreversibility ratio 


The method relies on the simultaneous application of principles 
of heat and mass transfer, fluid mechanics, and engineering ther- 
modynamics, in the pursuit of realistic models of processes, devices, 
and installations. 


2. A brief theory 


It is well known from the thermodynamics that heat transfer 
in the system results in thermodynamic irreversibility and it gen- 
erates the entropy. The minimization of entropy is a method of 
optimization of the thermodynamic imperfections and fluid flow 
irreversibilities [5]. Thus, investigation of entropy generation and 
minimization in the heating systems is important to include irre- 
versibilities into the thermodynamic analysis. Entropy is a measure 
of chaos and if there is entropy in a system, the quality of energy 
decreases. A system of high entropy is more chaotic or disordered 
than one of low entropy. For example, a field with papers scattered 
about has higher entropy than the field with the papers neatly piled 
[1]. A study of entropy generation in fundamental convective heat 
transfer was obtained by Bejan [6,7]. 

Entropy analysis involves the evaluation of the inputs, outputs, 
accumulations and creations of entropy for a considered system. It 
means that 


By adding mass transfer the equations [10], 
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In above equation, C defines the concentration. Eq. (3) clearly 
shows that entropy generation in a convection heat transfer 
occur mainly heat transfer and fluid friction. Thus, an important 
parameter can be defined as irreversibility ratio. The definition of 
irreversibility ratio may enhance the understanding of the irre- 


versibilities associated with the heat transfer and the fluid friction 
[5]. It is given as [6,7], 
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¢= a, (5) 
G, heattransfer 

Carrington and Sun [38] investigated the second law analysis 
of combined heat and mass transfer for both internal and exter- 
nal flows. They found an equation for the entropy generation. The 
dimensionless Bejan number (Be) [11] is used in entropy studies. It 
can be used as the alternative irreversibility distribution [11], 

2 
Be=— OD (6) 
(VT? + (VY? 

Be> 0.5 is the limit at which, heat transfer irreversibility dom- 
inates; Be «0.5 is the opposite limit at which, irreversibility is 
dominated by fluid friction effects; and Be = 0.5 is the case wherein, 
heat transfer irreversibility and fluid flow irreversibility are of equal 
importance. Here also, integrating Bejan number yields global Bejan 
number [12], 
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Entropy input + Entropy generation — Entropy output = Entropy accumulation change in the total entropy of the system)(1) 


In above formulation, the entropy creation is identical to negen- 
tropy consumption and proportional to exergy consumption [1]. 
Entropy balance in the rate form is given by 


ASsystem = (Sintet = Soie) T $ (2) 


The local volumetric rate of entropy generation, si (W/m? K), 
based on the second law of thermodynamics, is given for three 
dimensional flow as follows [8,9], 
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The first term on the right-hand side is due to heat transfer 
and the second term is due to viscous dissipation. It is well known 
that the viscous dissipation in the first law of thermodynamics is 
negligible in many engineering applications. 


Other fundamental studies on entropy generation in convec- 
tive heat transfer and isothermal convective mass transfer [13], 
combined heat and mass transfer [14], counter flow heat exchang- 
ers for gas-to-gas applications [15], fundamental convective heat 
transfer [16], design optimization of thermal system [17] and irre- 
versibilities in various duct [18]. Table 1 illustrates a summarize on 
fundamental studies of entropy generation. 


3. Studies for natural convection in porous media filled 
systems 


Natural convection in porous media has an important topic 
due to the day-by-day improvements of engineering in industry. 
There are many applications for this topic such as building materi- 
als, insulation technology, bioengineering, geothermal application, 
chemical technologies, filter design, etc. [19,20]. Entropy analysis is 
an effective method for studies on porous media to obtain system 
efficiency. Fundamentals of entropy generation in natural convec- 
tion are defined widely by Baytas and Baytas [21]. Mahmud and 
Fraser [22] investigated the nature of heat transfer and entropy 
generation for natural convection in a two-dimensional circular 
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Table 1 
Summary of fundamental studies on entropy generation in convective heat transfer. 
Author System Topic 
Rosen [1] Different engineering systems Second-law analysis 
Bejan [2] Different engineering systems Second-law analysis in heat transfer 
Bejan [6-9] Different engineering systems Fundamentals of exergy analysis, entropy generation minimization 


Arpaci and Selamet [27] 
Carrington and Sun [38] 


Energy systems 


external flows 
Bejan [15] 
gas-to-gas applications 
Design optimization 
Simulation 
Thermodynamics processes 


Tsatsaronis et al. [17] 
Naterer and Camberos 
Prigogine [63] 


Oztop [64] Semicircular duct 
Amani and Nobari [73] Curved pipes 
Chen et al. [67] Vertical channel 


Hernandez et al. [74] Fuel cell 


Combined heat and mass transfer in both internal and 


Heat exchanger design counter flow heat exchangers for 


Entropy efficiency 
Second-law analysis 


Irreversibility 


Exergy 

Entropy and the second law 
Irreversibilities 

Second law analysis 
Entropy generation 
Entropy generation 
Entropy generation analysis 


section enclosure filled with porous media. In their work, the Darcy 
momentum equation is used to model the porous media. Results 
are presented in terms of Nusselt number, entropy generation num- 
ber, and Bejan number. They found that for conduction regime, 
both average Nusselt number and the entropy generation number 
are independent of Rayleigh number variation. In convection- 
dominated regime, these parameters show an increasing tendency 
with increasing Rayleigh number. At high Rayleigh number, the 
near-wall magnitude of overall entropy generation rate is higher, 
but heat transfer irreversibility is higher at the center portion of the 
cavity. These are summarized in Fig. 1. By plotting Nusselt number 
and entropy generation number versus Rayleigh number. In their 
another study, the problem of entropy generation in a fluid sat- 
urated porous cavity for laminar magnetohydrodynamic natural 
convection heat transfer is analyzed by considering Darcy’s law for 
porous media. Magnetic force is assumed acting along the direc- 
tion of the gravity force. As boundary conditions of the cavity, two 
vertical opposite walls are kept at constant but different temper- 
atures and the remaining two walls are kept thermally insulated. 
For a range of Rayleigh number (Ra= 1-104) and Hartmann num- 
ber (Ha=0-10), heat transfer, overall entropy generation rate, and 
heat transfer irreversibility are presented in terms of dimension- 
less Nusselt number (Nu), entropy generation number (Ns), and 
Bejan number (Be), respectively. Fig. 2, in their paper, presents Nsay 
and Beay for Fig. 2(a) and (b), respectively. These are function of 
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Fig. 1. Average Nusselt and entropy generation numbers as a function of Rayleigh 
number [5]. 


Hartmann number at different Ra. Similar to the Nuay—Ha profile 
characteristics, Nsay is maximum at Ha =0, decreases with increas- 
ing Ha and approaches an asymptote at higher values of Hartmann 
number. Magnitudes of Beay are almost same for all Rayleigh num- 
bers considered here at Ha=0. Increases in the value of Ha have a 
tendency to slowdown the fluid movement inside the cavity, thus 
causing relative increases of heat transfer irreversibility (i.e., Beay). 

Baytas [11] performed a study to investigate the entropy gen- 
eration due to buoyancy forces in an inclined porous cavity. He 
observed that when Rayleigh number decreases, heat transfer irre- 
versibility begins to dominate the fluid friction irreversibility. The 
Bejan number is rapidly changed between 150° and 270°. The prob- 
lem of entropy generation for unsteady natural convection and 
radiation in a tilted open-ended channel, submitted to uniform 
hot walls temperature and filled with an isotropic fluid saturated 
porous medium is outlined by Slimi [23]. In his study, Darcy flow 
model was used and viscous dissipation effects are taken into 
account. His results showed that the global heat transfer rate 
depends slightly on the viscous dissipation term but is strongly sub- 
ordinate to the tilted angle. Moreover, the calculation of entropy 
generation maps is feasible and gives a measure of the degree 
of irreversibility of the convective flow. Recently, entropy gener- 
ation due to buoyancy-induced convection and conduction in a 
right angle trapezoidal enclosure filled with fluid saturated porous 
medium has been performed numerically by Varol et al. [19]. In 
their case, the left vertical solid wall of the trapezoidal enclosure 
has a finite thickness and conductivity. The governing Darcy and 
energy equations are solved numerically using a finite difference 
method. Entropy generation is calculated by using the obtained 
velocities and temperature distributions from the computer code. 
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Fig. 2. (a) Average entropy generation number as a function of Rayleigh number. 
(b) Average Bejan number as a function of Rayleigh number [16]. 
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Results are presented for the Bejan number, local and mean Nusselt 
numbers, streamlines, isotherms, iso-Bejan lines and entropy gen- 
eration contours. It is found that the most important parameters 
on heat transfer and fluid flow are thermal conductivity ratio and 
dimensionless thickness of the solid wall of the enclosure. Thus, 
these parameters also generate entropy for the whole system. It 
is also observed that increasing the Rayleigh number decreases 
the Bejan number; however, heat transfer is an increasing func- 
tion of Rayleigh number. Non-uniform heating effects are tested 
on entropy generation due to natural convection in a isosceles tri- 
angular enclosures at different positions by Varol et al. [24]. They 
found that both inclination angles and Rayleigh numbers make 
important effect on natural convection heat transfer, fluid flow and 
entropy generation. The highest entropy generation due to HTI and 
FFI and stream function are observed at ọ = 90°. Multiple cells were 
formed at this angle. Streamlines, isotherms and entropy contours 
are symmetric inside the enclosure for both g=0° and g=180°. 

In Zahmatkesh’s [12] study, a numerical analyze was per- 
formed on the importance of thermal boundary conditions of the 
heated/cooled walls in heat transfer and entropy generation char- 
acteristics inside a porous enclosure, heated from below. Both 
the heating and the cooling are carried out uniformly and non- 
uniformly and the results are compared. He presents results in 
terms of streamlines, isothermal lines, iso-entropy generation lines, 
and iso-Bejan lines. Additionally, variations of average Nusselt 
number, global entropy generation rate, and global Bejan number 
are analyzed over a wide range of Darcy-modified Rayleigh num- 
ber (10 < Ra < 1000). Inspection of the results indicates that thermal 
boundary conditions are of profound influences on the induced 
flow as well as heat transfer characteristics and possess prominent 
consequences on entropy generation rates. It is demonstrated that, 
the optimum case with respect to heat transfer as well as entropy 
generation could be achieved by non-uniform heating. 

In another study of Baytas [25], the entropy generation in a sat- 
urated porous cavity for laminar natural convection is investigated 
by using the non-Darcy flow model and a thermal nonequilibrium 
model of the heat transfer between the fluid and the solid phases. 
His results demonstrated that the entropy generation due to vis- 
cous drag is very important for the non-Darcy flow model, and the 
high temperature difference between solid and fluid phases shows 
the importance of the thermal nonequilibrium model for porous 
media. He suggested that entropy generation for the thermal 
nonequilibrium and non-Darcy flow models should be considered 
for engineering design and thermodynamic optimization. Mahmud 
and Fraser [26] investigated the nature of entropy generation for 
natural convection in a square enclosure filled with porous media 
vibrating sinusoidally perpendicular to the applied temperature 
gradient in a zero-gravity field. They indicated that Gravity oscil- 
lation introduces a true periodic behavior to the Nusselt number, 
Bejan number, and entropy generation rate. The periodic response 
of these three parameters is synchronized with the forced acceler- 
ation, namely, having the same period as the forced acceleration. 
Table 2 gives a summary of the various theoretical investigations 
for entropy generation in porous media filled systems. 


4. Studies for natural convection viscous fluid filled 
medium 


In industrial application, natural fluids as air or water are used 
mostly. Thus, in the analyses, thermal properties of these fluids 
are taken into account to the analysis. Arpaci and Selamet [27] 
investigated the internal energy generation in case of buoyant tur- 
bulent flows. They developed a thermal model based on microscales 
flows for different parameters. Dagtekin et al. [28] performed a 
numerical study to investigate the second law analysis due to 
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Fig. 3. Evolution of the maximum of the dimensionless total entropy genera- 
tion versus Rayleigh number and distribution irreversibility ratio g=10-4, 10-3, 
3.3 x 10-3 and 10-? [30]. 


buoyancy-induced flow in a I’-shaped enclosure. They investigated 
effects of geometrical ratio and Rayleigh number on entropy gen- 
eration. 

Influence of thermodiffusion on the birth of free convection 
in Rayleigh-Benard configuration was investigated by Platten and 
Chavepeyer [29]. They summarize the influence of thermodiffusion 
on free convection together with the way to use this mutual influ- 
ence in order to experimentally deduce the value of the Soret effect. 
Magherbi et al. [30] studied the time-dependent natural convec- 
tion and second law analysis of thermodynamics in a differentially 
heated enclosure using finite volume technique. They indicated 
that the total entropy generation has a maximum value at the 
onset of the transient state, which increases with the Rayleigh num- 
ber and the irreversibility distribution ratio. They also found that 
entropy generation asymptotically tends towards a constant value 
at low Rayleigh numbers, whereas an oscillation of the entropy gen- 
eration was observed for higher Rayleigh numbers, before reaching 
the steady state. Graphs of the maximum of the total entropy 
generation Max (Sra) as function of Rayleigh number at differ- 
ent irreversibility distribution ratio g are shown in Fig. 3. At low 
irreversibility distribution ratios (as low as y < 1074) Max (Sra) 
takes on small values even at high Rayleigh numbers. However, 
for y < 10-3, Max (r a) increases rapidly. Similar observations can 
be made about the evolution of the total entropy generation (Sra) 
as a function of the Rayleigh number during the steady state as seen 
in Fig. 4. 

Entropy generation in rectangular cavities with the same area 
but different aspect ratios is numerically investigated by Ilis et al. 
[31]. They used finite difference method and the energy and vortic- 
ity equations are solved line by line by employing the ADI method, 
whereas the stream function equation is solved point by point. It 
is found that for a cavity with high value of Rayleigh number (i.e., 
Ra=10°), the total entropy generation due to fluid friction and total 
entropy generation number increase with increasing aspect ratio, 
attain a maximum and then decrease. For all cavities with different 
aspect ratio, the total entropy generation increases with increasing 
Rayleigh number, though the rate of increase is different as seen 
in Fig. 5. For Ra=100, the total entropy generation increases from 
the minimum aspect ratio (AR=1) to the maximum aspect ratio 
(AR=16) with the same order. For Ra=10° and ¢=10~, the total 
entropy generation for the AR = 16 cavity is less than that for AR=9 
due to the peak point of the fluid friction total entropy generation. 

Mourad et al. [10] report a numerical determination of the 
entropy generation in doubly diffusive convection on 2D approx- 
imation in a square inclined cavity, filled with a fluid (assumed 
to be a perfect gas mixture). Their results showed that the total 
entropy generation increases with the thermal Grashof number 
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Table 2 
Summary of investigations in porous media filled enclosures. 


Author Study 


Method 


Results and remarks 


Mahmud and Fraser [16] Magnetohydrodynamic free 
convection in a porous 
enclosure 

Free convection inside a 
two-dimensional circular 


porous enclosure 


Mahmud and Fraser [22] 


Zahmatkesh [12] Free convection inside a 
circular porous enclosure with 


different boundary conditions 


Baytas [25] Porous square enclosure filled 
with a heat-generating solid 
phase 

Baytas [11] Tilted saturated porous cavity 


Control volume based 
finite-volume method 


Control volume based 
finite-volume method was 
used. A non-staggered and 
non-uniform grid system is 
used with a higher mesh 
density near the walls 
Control volume method of 
Patankar [62] 


Non-equilibrium porous model 
Finite volume finite difference 
method 

Finite difference control 
method of Patankar [62] 

Finite difference method 


In the absence of magnetic force, entropy generation rate 
is relatively higher in magnitude near two vertical walls 


At high Rayleigh number, the near-wall magnitude of 
overall entropy generation rate is higher, but heat transfer 
irreversibility is higher at the center portion of the cavity 


Entropy generation rate is likely to be the highest for 
uniform heating/cooling and the lowest for non-uniform 
heating 

The entropy generation due to viscous drag is very 
important for non-Darcy flow model as well as high 
temperature difference 

Inclination angle of cavity makes important effect on 
degree of irreversibility 


Varol et al. [19] Right-angle trapezoidal 


enclosure 
Slimi [23] Tilted saturated porous 
channel 
Mahmud and Fraser [49] Square porous enclosure with 
vibration volume method 


Control volume based finite 


Inclination angle of the inclined wall of trapezoidal 
geometry affects the entropy generation 

Entropy generation maps is feasible and gives a measure of 
the degree of irreversibility of the convective flow 

At the lower extreme of the gravity oscillation Nua, and 
Nsay is minimum and Bey is maximum 


and the buoyancy ratio for moderate Lewis numbers. Locally, the 
irreversibility due to heat and mass transfer are nearly identical and 
are localized in the bottom and the top of the heated and the cooled 
walls respectively. Entropy generation due to natural convection 
has been calculated for three radii and a wide range of Rayleigh 
numbers for an isothermal cylinder by Abu-Hijleh et al. [32]. They 
indicated that, for air, entropy generation was predominantly due 
to conduction while the contribution of viscous dissipation was 
negligible. As the Rayleigh number increased, total entropy genera- 
tion increased. As the radius of the cylinder increased, total entropy 
generation decreased. 
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Fig.4. Variation of the dimensionless total entropy generation in steady state versus 
Rayleigh number and distribution irreversibility ratio y= 10-4, 10-3, 3.3 x 10-3, and 
10-2 [30]. 


The optimization in an inclined square enclosure for minimum 
entropy generation was analyzed by Baytas [33]. Based on that 
study, the local heat transfer irreversibility and the local fluid 
friction irreversibility change by the inclination angle and the min- 
imum entropy generation depends considerably on the inclination. 
Numerical prediction of local and total entropy generation rates for 
natural convection of air in a vertical channel symmetrically heated 
with a uniform heat flux was studied by Andreozzi et al. [34]. They 
found that global entropy generation increases with both aspect 
ratio and Rayleigh number increase. These values have been corre- 
lated as function of Ra and L/b; the correlation obtained by means 
of the least square method is 


log S% y = -13.2 + 1.91 log Ra + 0.045 (5) (8) 
In above equations, 
r? = 0.999 and 5< ; <20 x 10? < Ra < 10°. (9) 


Erbay et al. [35] analyzed the thermodynamic irreversibilities 
due to natural convection in a square enclosure. Erbay et al. [36,37] 
investigated entropy generation with a single discrete heat source 
on the vertical wall of a cavity. Their work focused on the transient 
development of the flow and temperature fields and instantaneous 
entropy generation distributions during the transient stages. 

Shuja et al. [5] made a numerical work on natural convection in 
a square cavity with a heat generating body by including entropy 
consideration. They solved both convection and conduction (conju- 
gate) mode of heat transfer to simulate cooling of micro-electronic 
equipments. They compared the effects of working fluid as air and 
water and found that the flow developed in the cavity due to nat- 
ural convection has almost the same mode of velocity fields for 
air and water. Cooling performance of heated body is related with 
location of the body. The entropy generation due to heat transfer is 
relatively lower for the solid body located in left bottom and right 
top corners of the cavity. It attains high values in the case of fluid 
friction. It is an interesting result that the entropy attains higher 
values in the case of air as compared to the case for water. It means 
that the heated solid body losses more heat in air than in water 
for the indicated simulation conditions. Table 3 illustrates the total 
entropy generation in the cavity due to HTI and FFI. In The case of 
air total entropy generation due to heat transfer and fluid friction 
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Table 3 
Total entropy generated in the cavity for various cases [45]. 


Working Fluid Solid position Entropy generated (J/K) due to 
Heat transfer Fluid friction Total Irreversibility ratio (p) 

Air Bottom 2.379 x 10-9 6.377 x 107!° 3.017 x 107° 0.2680 

Center 3.228 x 10-9 1.286 x 10-10 3.357 x 10-9 0.0399 

Top 2.419 x 10-9 2.797 x 10-10 2.698 x 10-9 0.1160 
Water Bottom 3.238 x 10-1! 1.170 x 10-1° 1.494 x 10-19 3.6100 

Center 5.884 x 107-1! 2.020 x 10-1! 7.905 x 107! 0.3430 

Top 2.560 x 10-1! 5.024 x 10-1! 7.584 x 10-1! 1.9600 


attains higher values as compared to that for water. In addition, the 
total entropy generation due to heat transfer has maximum value 
when the solid body is at the center of the cavity while the total 
entropy generation due to fluid friction has a minimum value at 
this location. This may indicate that fluid motion generated close 
to the solid body is equally distributed in the vicinity of all sur- 
faces of the solid body, which results in a relatively low velocity 
flow field. On the other hand, the increase in total entropy due to 
heat transfer is not only because of the convection effect, but the 


0 - ; 
10° 10' 10° 10° 10° 
Ra 


Fig. 5. The changes of total entropy generation with Rayleigh number for cavities 
with different aspect ratios (a) ®=10~4 and (b) =10-? [31]. 


conduction effect as well. From the table, the irreversibility ratio is 
minimum when the solid body is located at the center of the cavity. 

Varol et al. [39] studied the forecasting of entropy generation 
of laminar natural convection in a partially cooled square cross- 
sectional room has been performed using support vector machines 
(SVM). The two-dimensional room was modeled as floor heating 
story with a window. Forecasting of entropy generation due to fluid 
friction irreversibility (FFI) and heat transfer irreversibility (HTI) 
were made with known values for unknown parameters using SVM. 
Thus, calculation time was extremely reduced and values were 
obtained even for non-convergence cases. 

Recently, Mukhopadhyay [40] made a numerical work to ana- 
lyze entropy generation due to natural convection in square 
enclosures with multiple discrete heat sources. He examixed the 
heaters of equal length and strength, the effects of Rayleigh num- 
ber and heater position on flow and temperature fields and local 
entropy generation. He indicated that for heaters of equal length 
and strength ratios, placement of the heaters close to the mid- 
plane of the cavity is most undesirable as both peak temperature 
and entropy generation are maximized. On the other hand, placing 
the heaters very close to the walls or at a dimensionless spacing of 
about 0.5 is most desirable as both these quantities are minimized. 

Yang et al. [41] performed the second law analysis of thermody- 
namics on the laminar film condensation of pure saturated vapor 
flowing in the direction of gravity on an ellipsoid with variable wall 
temperature. They found that entropy generation increases with 
ellipticity. Furthermore, the irreversibility due to finite tempera- 
ture difference heat transfer dominates over that due to condensate 
film flow friction and the local entropy generation rate decreases 
with increasing ellipticity in the upper half of ellipsoid. Abu-Hijleh 
et al. [42] made a numerical prediction of entropy generation due 
to natural convection from a horizontal cylinder. Fig. 5 shows that 
an increase in Rap, i.e., a higher temperature difference between 
the cylinder and surroundings, results in an increase in Sgen. The 
main conclusion from Fig. 5 is that the use of larger cylinders 
results in significantly lower entropy generation for the same value 
of Rap due to a lower temperature difference. A larger cylinder 
has a lower tangential temperature gradient, which reduces the 
contribution from the second component in the conduction part 
of Eq. (8). Increased cylinder size results in an increase in vis- 
cous dissipation in Eq. (8). It should, however, be noted that the 
second term of Eq. (8) was negligible even at the largest ro for 
the greatest Rap. Thus, there was no appreciable increase in Sgen 
at large ro due to viscous dissipation. Analysis of local entropy- 
generation profiles showed that the location of maximum entropy 
generation depends on the size of the cylinder. For rọ =10-3 m, 
maximum entropy generation occurred in the plume region, at 
the top of the cylinder. For rg = 10-! m, maximum entropy genera- 
tion occurred in the inflow region, at the bottom of the cylinder. 
For 107! >ro > 1073, the contributions from both the inflow and 
plum regions were of the same order. We may manipulate the flow 
in regions of low local entropy generation with minimal adverse 
effects on the total irreversibility of the component. Jery et al. [43] 
performed a study on influence of an external oriented magnetic 


Table 4 


H.F. Oztop, K. Al-Salem / Renewable and Sustainable Energy Reviews 16 (2012) 911-920 917 


Summary of the calculation of entropy due to natural convection in different enclosures. 


Authors 


Geometry 


Method 


Results and remarks 


Abu-Hijleh and Heilen [46] 


Abu-Hijleh et al. [32] 


Abu-Hijleh et al. [42] 


Andreozzi et al. [34] 


Baytas [33] 


Mukhopadhyay [40] 


Famouri and Hooman [52] 


ilis et al. [31] 


Magherbi et al. [30] 


Mourad et al. [10] 


Erbay et al. [37] 


Shuja et al. [5] 

Mahmud and Sadrul Islam [53] 
Varol et al. [54] 

Varol et al. [55] 


Mahmud and Fraser [56] 


Bassam [58] 
Yilbas [57] 


Mahmud et al. [59] 
Oliveski et al. [69] 
Kaluri and Basak [68] 
Saleem et al. [70] 


Kaluri and Basak [71] 


Basak et al. [72] 


Heated rotating cylinder 


Heated isothermal cylinder 
in oil 


Natural convection from a 
horizontal cylinder 


Vertical channel 
symmetrically heated at 
uniform heat flux 


Natural convection in 
inclined cavity 

Natural convection in 
square enclosures with 2 
multiple discrete heat 
sources 

Natural convection by 
heated partitions ina 
cavity 


Rectangular cavity with 
differentially heated 
vertical walls 


Entropy generation at the 
onset of natural convection 
in a differentially heated 
enclosure 

Numerical determination 
of the entropy generation 
in doubly diffusive 
convection on 2D 
approximation in a square 
inclined cavity 

An analysis of the entropy 
generation in a square 
enclosure with partial 
heating 

Square cavity with heat 
generating body 

Wavy enclosure 


Partially heated isosceles 
triangular enclosures 
Conjugate natural 
convection in enclosures 
Mixed 
convection-radiation 
interaction in a vertical 
channel 

Cylinder 

Cavity heated from bottom 


Vertical channel 


Enclosures 
Differentially heated 
cavities 

Marangoni convection 


Discretely heated square 
cavity 


Wide variety of thermal 
boundary conditions 


Elliptic system of PDEs are 
solved by explicit hybrid 
scheme 

Elliptic system of PDEs are 
solved by explicit hybrid 
scheme 

Equations are solved in 
streamfunction-vorticity 
formulation 

Numerical solutions of the 
elliptic momentum and 
energy equations are 
carried out with the 
stream-function-vorticity 
method 

Finite volume technique 


The governing equations 
were discretized using 
SIMPLER algorithm of 
Patankar [62] 

FORTRAN code and the 
commercially available 
CFD-ACE software. 


Finite difference technique 


Of the control volume 
finite-element method 
(CVFEM) 

Control volume 


finite-element method 
(CVFEM) 


Finite volume method 
(FVM) coupled with power 
law scheme for the 
convective terms 

Finite volume method 
Finite difference method 
Finite difference method 
Analytical study with 
reasonable simplifications 
Streamfunction-vorticity 
Finite volume method 
Analytical solution 
General methods 

Finite element method 
Finite difference method 
Finite element method 


Galerkin finite element 
method 


The total entropy generation increases with the increase of 
both Reynolds number and buoyancy parameter but 
decreases as the cylinder radius increases 

Including the viscous dissipation term had minimal effect 
on Nusselt number and total entropy generation 


As the radius of the cylinder increased, total entropy 
generation decreased. 


Global entropy generation increases with both aspect ratio 
and Rayleigh number increase 


An optimization has been made to obtain minimum 
entropy generation 

The dominant contribution to entropy generation comes 
from heat transfer irreversibilities, with fluid friction 
accounting for only a small fraction of the total entropy 
generation 

While fluid friction term has nearly no contribution to 
entropy production, the heat transfer irreversibility 
increases monotonically with the Nusselt number and the 
dimensionless temperature difference 

For a cavity with high value of Rayleigh number (i.e., 
Ra=105), the total entropy generation due to fluid friction 
and total entropy generation number increase with 
increasing aspect ratio, attain a maximum and then 
decrease. 

The total entropy generation has a maximum value at the 
onset of the transient state, which increases with the 
Rayleigh number and the irreversibility distribution ratio 


Total entropy generation increases with the thermal 
Grashof number and the buoyancy ratio for moderate 
Lewis numbers 


The irreversibilities are dominant due to heat transfer 
whereas fluid friction irreversibilities have been found 
negligible as it is expected for the natural convection 


The entropy generation due to heat transfer is relatively 
lower 

Wave ratio is an important parameter on entropy 
production 

Length of partial heater the most effective parameter on 
entropy generation 

Thermal conductivity ratios and thickness of solid wall 
become effective on entropy production 

Irreversibility due to the heat transfer and/or fluid friction 
become zero at idle points 


Entropy increases with Ra number 

Temperature difference play important role on entropy 
generation 

For positive value of heat generation/absorption 
parameter, entropy generation rate is higher than the 
negative value of same magnitude 

Entropy generation is applied for natural convection 
Prandtl number plays important role on entropy 
generation 

Entropy generation rate increases with the increase in the 
Marangoni number 

High thermal mixing may not be the optimal situation for 
achieving uniform temperature distribution based on 
entropy production 

Thermal boundary conditions are affected on thermal 
mixing, heat transfer and entropy generation 
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field on entropy generation in natural convection for air and liq- 
uid. They studied the steady-unsteady states mass, momentum and 
entropy equation for different inclination angle of the magnetic 
field and Pradntl number. At local level and for relatively higher 
thermal Grashof number (Grr = 10°), entropy generation distribu- 
tion is strongly dependent on magnetic field direction, magnitude 
of irreversibility lines increases up to 30°, then gradually decreases. 
No entropy is generated in the cavity center. Table 4 summarizes 
the various investigations for entropy generation in viscous fluid 
filled systems. 


5. Entropy generation due to mixed convection 


Mixed convection heat transfer occurs for both natural and 
forced convection present in a system. In this case, Richardson 
number (Gr/Re?) becomes effective. Entropy places important role 
in this type of heat transfer regime due to heat transfer and fluid 
friction. 

Narusawa [44] studied the second-law analysis in a rectangular 
duct for mixed convection regime. The simulation of mixed con- 
vection in a square cavity with protruding body having different 
aspect ratios was performed by Yilbas et al. [45]. Abu-Hijleh and 
Heilen [46] studied the entropy generation due to laminar mixed 
convection from an isothermal rotating cylinder for three cylinder 
radii and covered wide ranges of Reynolds number and buoyancy 
parameter. They observed that entropy generation increased as 
the Reynolds number and buoyancy parameter increased and it 
decreased as the cylinder radius increased. For the same com- 
bination of Reynolds number and buoyancy parameter, entropy 
generation was mainly due to thermal effects at small cylinder radii 
and due to viscous effects at large cylinder radii. Abbassi et al. [47] 
analyzed the entropy generation in Poiseuille—Benard channel flow. 
He found that the maximum entropy generation is localized at areas 
where heat exchanged between the walls and the flow is maximum. 
There is no significant entropy generation is seen in the main flow. 
In the similar manner, Nourullahi et al. [48] interested by the issue 
of entropy generation and Nusselt number in Poiseuille-Benard 
channel flow are analyzed by solving numerically Navier-Stokes 
and energy equations with the use of the classic Boussinesq in 
compressible approximation. Their result showed that the Nus- 
selt number changes very slightly and it is almost constant for 
low values of inclination angle. It decreases for higher values of 
inclination angle. The entropy generation due to heat transfer is 
localized at areas where heat exchanged between the walls and 
the flow has a maximum value, while the entropy generation due 
to fluid friction is maximum at areas where the velocity gradients 
are maximum such as vortex centers. In other words, the inclina- 
tion angle is a good parameter to control entropy generation in the 
presence of buoyancy forces. Mahmud and Fraser [49] focused on 
magnetohydrodynamic mixed convection through a vertical chan- 
nel packed with fluid saturated porous substances by applying first 
and second laws of thermodynamics to analyze the problem. They 
solved governing equations analytically under reseasonable sim- 
plifications for different parameters such as Hartmann number, 
Plank number, Richardson number and group parameter (Br//7). 
They found that the entropy generation number is characterized 
by a concave shaped profile and is symmetric about the chan- 
nel centerline for a symmetrical temperature boundary condition. 
Radiation and mixed convection parameters have a more domi- 
nating influence on entropy generation rate than porous-magnet 
and group parameters. Expressions for the idle points of entropy 
generation are derived and their location(s) are determined. Irre- 
versibility due to the heat transfer and/or fluid friction becomes 
zero at these idle points. Shohel and Roydon [50] did a similar work 
on mixed convection-radiation interaction in a vertical channel by 


calculating entropy generation. Another application has been per- 
formed by Cheng et al. [51] on entropy production in mixed 
convective flow in a vertical channel with transverse fin arrays. 


6. Entropy generation of natural convection in nanofluid 
filled enclosures 


Addition of nanoparticles (Al203, Cu, Ag, Au, CuO, TiO}, etc.) of 
any material into base fluid such as water, glycerol, etc. is called 
as nanofluid. In this way, heat transfer of fluid is enhanced due to 
higher thermal conductivity values of nanoparticles. In the same 
time, viscosity of fluid is changing and it affects the heat trans- 
fer enhancement. Studies on entropy generation of heat transfer 
in nanofluid filled energy systems are extremely limited. In this 
context, Shahi et al. [65] performed a numerical study on entropy 
generation due to natural convection of a nanofluid that consists 
of water and Cu in a cavity with a protruded heat source. They 
used finite volume method and solved the equations for different 
Rayleigh numbers, solid concentration and heat source location. 
Their results showed that the maximum value of Nusselt number 
and minimum entropy generation are obtained when heat source 
mountains in the bottom horizontal wall. Singh et al. [66] present a 
theoretical investigation of the entropy generation analysis due to 
flow and heat transfer in nanofluids. They considered the most com- 
mon alumina—water nanofluids as the model fluid. They observed 
that at lower tube diameter, flow friction irreversibility is more sig- 
nificant and at higher tube diameter thermal irreversibility is more. 
For both laminar and turbulent flow, there is an optimum diameter 
at which the entropy generation rate is the minimum for a given 
nanofluid. 


7. Conclusions 


The present review is a comprehensive outlook on the research 
progress made on entropy generation in natural convection and 
mixed convection in ducts or channels. The important result can 
be drawn from the reviewed literature is that the main reasons of 
entropy generation are heat transfer and fluid friction in a thermal 
system with buoyancy force. From the observed results it is clearly 
seen, that analysis of entropy generation is crucial for all thermal 
systems to reduce energy and control of heat transfer and fluid 
flow. The most important parameter is Rayleigh number on entropy 
generation in natural convection. The entropy generation increases 
with increasing of Rayleigh number due to increasing of heat trans- 
fer and fluid friction. Entropy analysis indicates a local address 
that where energy losses high in the physical model or systems. 
Darcy number is also another effective parameter for entropy gen- 
eration in porous media filled systems. The entropy also increases 
with increasing of Darcy number. In the case of mixed convection, 
the entropy generation decreases with increasing of Richardson 
number due to decreasing of heat transfer. Hartmann number 
is another effective parameter in magnetohydrodynamics studies 
with electrically conductive fluids. In this case, increasing of Hart- 
mann number decreases the flow velocity. However, the entropy is 
reduced. Application direction of magnetic force is also important 
to minimize the entropy generation. Entropy generation decreases 
with decreasing of irreversibility ratio. Inclination angle makes 
important effect in natural convection studies due to heat transport 
way changes by gravity. In these kinds of studies, entropy genera- 
tion is also affected. The geometrical shape and aspect ratio of the 
cavity change the entropy generation even at constant Rayleigh 
number. And type of the fluid with Prandtl number places impor- 
tant role. Recently, the second law thermodynamics was applied 
to industrial problems to minimize the entropy generation [60,61]. 
There are very few correlations between effective parameters and 
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entropy generation. Researcher can focus on that subject for future 
studies. The trend in the literature goes to extension of natural and 
mixed convection studies by adding entropy generation. 
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